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Due  to  mismatch  of  the  coefficients  of  thermal  expansion  of  various  layers  in  the  posi¬ 
tive/electrolyte/negative  (PEN)  structures  of  solid  oxide  fuel  cells  (SOFC),  thermal  stresses  and  warpage 
on  the  PEN  are  unavoidable  due  to  the  temperature  changes  from  the  stress-free  sintering  temperature 
to  room  temperature  during  the  PEN  manufacturing  process.  In  the  meantime,  additional  mechanical 
stresses  will  also  be  created  by  mechanical  flattening  during  the  stack  assembly  process.  In  order  to 
ensure  the  structural  integrity  of  the  cell  and  stack  of  SOFC,  it  is  necessary  to  develop  failure  criteria  for 
SOFC  PEN  structures  based  on  the  initial  flaws  occurred  during  cell  sintering  and  stack  assembly.  In  this 
paper,  the  global  relationship  between  the  critical  energy  release  rate  and  critical  curvature  and  maximum 
displacement  of  the  warped  cells  caused  by  the  temperature  changes  as  well  as  mechanical  flattening 
process  is  established  so  that  possible  failure  of  SOFC  PEN  structures  may  be  predicted  deterministically 
by  the  measurement  of  the  curvature  and  displacement  of  the  warped  cells. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cells  are  high-efficiency  energy  conversion  devices  that  are 
environmentally  friendly  with  little  or  no  toxic  emission.  With  high 
efficiency  and  potential  applications  in  stationary  power  generation 
as  well  as  auxiliary  power  units,  solid  oxide  fuel  cell  (SOFC)  contin¬ 
ues  to  show  great  promise  as  a  future  power  source.  Among  various 
SOFC  designs,  anode-supported  planar  cells  seem  to  provide  the 
best  performance  at  the  reasonable  cost  [1,2].  In  anode-supported 
cells,  a  porous  anode  provides  the  main  structural  support  for  the 
thin  (~10  pan)  yttria-stabilized  zirconia  (YSZ)  electrolyte  layer.  The 
cells  must  withstand  a  significant  number  of  thermal  cycles  with 
the  build-up  of  internal  stresses,  which  may  lead  to  the  cracking  of 
the  electrolyte  or  delamination  of  the  electrolyte-cathode/anode 
interface  [3]. 

Due  to  mismatch  of  the  coefficient  of  thermal  expansion  (CTE) 
of  various  layers  in  the  positive/electrolyte/negative  (PEN)  struc¬ 
tures  of  SOFC,  thermal  stresses  on  the  PEN  are  unavoidable  due  to 
the  temperature  changes  from  the  stress-free  sintering  tempera¬ 
ture  to  room  temperature  during  the  PEN  manufacturing  process. 
The  stresses  in  various  layers,  such  as  anode,  electrolyte,  and  cath¬ 
ode,  may  cause  the  initiation  and  propagation  of  cracks.  After  the 
sintering  processes  of  PEN  structure,  cracks  were  observed  in  PEN 
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structures  of  SOFC.  The  cracks  were  generated  between  the  elec¬ 
trolyte  (YSZ)  and  fuel  electrode  (NiO-YSZ)  in  PEN  structure  of 
fuel  cells  [4,5].  Such  thermal  internal  stresses  will  also  cause  cell 
warpage  with  a  certain  curvature.  Although  various  techniques 
have  been  developed  for  sintering  the  cells,  some  degree  of  cell 
warpage  remains  after  sintering.  During  stack  assembly,  the  cells 
with  warpage  will  be  flattened,  which  will  contribute  to  the  internal 
stresses  and  the  initiation  and  propagation  of  cracking  and  delam¬ 
ination.  In  [6],  it  was  reported  that  processing-induced  residual 
stresses  in  the  component  layers  of  SOFCs  can  lead  to  fracture  or 
to  cell  curvature  which  impedes  stack  assembly.  Generally,  it  has 
been  observed  that  there  are  two  basic  fracture  mechanisms  occur¬ 
ring  in  SOFC  stacks:  delamination  along  the  interface  and  transverse 
cracking  across  the  layers  [7].  Given  the  severe  thermal  cycles  that 
SOFCs  must  undergo,  it  is  essential  to  develop  suitable  materials  and 
design  stacks  so  that  they  are  able  to  sustain  thermo-mechanical 
loads  without  cracking. 

To  ensure  the  structural  integrity  of  the  SOFC  stack  without 
cracking  and  fracture,  many  efforts  have  been  paid  on  investigat¬ 
ing  the  fracture  behavior  of  SOFCs  and  its  PEN  with  the  various 
cracks  and  flaws  by  numerical  analyses  [8-10].  In  the  numerical 
fracture  analysis,  pre-crack  has  to  be  formed  in  the  various  inter¬ 
faces  of  PEN  structures  and  inside  various  layers,  and  fine  mesh 
has  to  be  used  at  the  crack  tips  to  capture  the  stress  field  of  the 
crack  tips  in  order  to  obtain  the  accurate  energy  release  rate.  3D- 
curved  crack  in  PEN  structures  was  investigated  by  means  of  the 
volume  integrals  model  in  [  10].  It  may  be  noted  that  the  SOFC  PEN  is 
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Nomenclature 

E  Young’s  modulus 

G  energy  release  rate 

Gc  interfacial  fracture  toughness 

H  anode  thickness 

L  nominal  cell  size 

M  flattening  moment 

N  resultant  forces  acting  on  all  the  layers 

M  resultant  moments  acting  on  all  the  layers 

Q.  elastic  constant  matrix 

Q  distributed  transverse  loading 

R  curvature  radius 

AT  temperature  difference 

u  layer  displacement  in  the  x-direction 

v  layer  displacement  in  the  y-direction 

w  layer  displacement  in  the  z-direction 

Uq  displacement  of  the  midplane  in  the  x-direction 

v0  displacement  of  the  midplane  in  the  y-direction 

w0  displacement  of  the  midplane  in  the  z-direction 

W  maximum  allowable  warpage 

Greek  symbols 

a  coefficient  of  thermal  expansion 

P  Dundurs  parameter 

k  warped  cell  curvature 

v  Poisson’s  ratio 


characterized  by  the  multilayer  thin  film  structure.  A  lot  of  efforts 
have  been  contributed  in  the  fracture  mechanism  of  the  thin  film 
structures.  Interfacial  fracture  toughness  of  the  thin  film  structures 
was  experimentally  and  theoretically  studied  in  [11-16].  The  chan¬ 
neling  cracks  of  the  thin  film  and  the  damage  of  its  propagation 
on  the  substrates  were  investigated  deeply  for  elastic  substrate  and 
elasto-plastic  substrates  in  [17-21  ].  The  detailed  review  is  referred 
to  [22].  In  these  works,  energy  release  rate  and  stress  intensity  fac¬ 
tors  are  widely  used  as  criteria  of  fracture  failure,  based  on  the 
local  stress  state  at  the  crack  tip  obtained  by  the  theoretical  and 
numerical  analyses. 

While  being  an  effective  method  for  the  reliability  design 
of  SOFC  PEN  structures,  the  fracture  toughness-based  approach 
described  above  will  have  difficulties  in  predicting  failure  induced 
during  the  sintering  and  mechanical  flattening  processes  for  the 
SOFC  PEN  materials.  This  is  because  the  energy  release  rate  may  not 
be  measured  directly  for  the  PEN  structures  of  SOFCs  during  sinter¬ 
ing  and  flattening  processes.  To  ensure  the  structural  integrity  and 
reliability  of  the  SOFC  stack,  it  is  necessary  to  develop  a  failure  cri¬ 
terion  for  the  initial  failures  that  occur  during  cell-stack  assembly. 

In  this  paper,  the  global  failure  criteria  for  PEN  fracture  are  estab¬ 
lished  based  on  local  fracture  mechanics  formulation  for  cracks  at 
various  locations  in  the  PEN  that  could  be  natural  result  on  var¬ 
ious  interfaces  and  in  cell  materials  due  to  the  porous  nature  of 
the  anode  and  cathode.  The  critical  energy  release  rate  is  related  to 
critical  curvature  and  maximum  deformation  of  the  warped  cells 
during  cell  fabrication;  therefore  the  allowable  curvature  of  the 
warped  cells  may  be  calculated,  and  possible  structural  failures  of 
the  PEN  structure  can  be  predicted  prior  to  stack  assembly.  Specifi¬ 
cally,  the  global  failure  criteria  are  formulated  based  on  the  warpage 
of  the  tri-layer  PEN  structure,  i.e.,  the  cell  after  sintering  process. 
The  warpage  of  each  cell  is  measured  before  stack  assembly.  Then, 
using  the  global  failure  criteria,  the  engineers  can  predict  whether  a 
cell  can  survive  the  stacking  assembly  process.  Clearly,  such  global 
failure  criteria  can  also  be  used  for  designing  the  sintering  pro¬ 
cesses  to  avoid  excess  warpage  of  the  cells  by  controlling  the  pore 
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Fig.  1.  A  functionally  graded  multiple  layer  laminate. 


size  and  porosity  using  the  sintering  temperature  and  time  as  well 
organic  compounds  in  the  PEN  structures.  The  purposes  of  develop¬ 
ing  such  global  failure  criteria  are  to  aid  the  initial  design,  materials 
selection,  and  optimized  assembly  parameters  of  SOFC  PENs. 


2.  Basic  governing  equations  of  thermal-mechanical 
analysis  of  multi-layers  and  global  fracture  criterion 

2  A.  Basic  governing  equations  of  thermal-mechanical  analysis 
for  functionally  graded  multi-layers 


Consider  a  multilayered  laminate  of  Lx  x  Ly,  as  shown  in  Fig.  1. 
The  cross-section  of  the  laminate  is  schematically  shown  in  Fig.  2. 
The  thermo-mechanical  properties  of  each  layer  may  vary  layer  by 
layer. 

Let  u,  v  and  w  be  the  displacements  in  the  x-direction, y-direction 
and  z-direction,  respectively.  Assuming  classical  plate  bending  the¬ 
ory  holds  and  that  plane  cross-section  remains  plane,  the  strain  in 
the  laminate  can  be  written  as  [23]: 


du  duo  d2w0 

Sx  =  =  ~z~d^’ 

dv  dv0  92w0 

Sy  =  dy  =  ~dy  ~Z~dp~’ 

du  dv  du0  dv0 

Vxy=ty  +  te  =  W  +  ^~ 


2  z 


d2w0 
dxdy  ' 


(1) 

(2) 

0) 


where  u0,  v0  and  w0  are  the  displacements  of  the  middle  plane  z  =  0. 
The  preceding  strain-displacement  relationship  can  be  written  in 
terms  of  the  midplane  strains  and  the  plate  curvature  as  follows: 


£  =  £°  —  ZK 


(4) 


Fig.  2.  Cross-section  of  a  multilayer  laminate. 
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where 


s°  =  [£X>  ey’  S°]T  = 


duo  dvo  3uo  dvo 
dx  ’  dy  ’  3y  +  dx 


k  =  [/cx,  /cy 


^xyj  — 


92w0  92w0  232w0 

3x2  ’  9y2  ’  9x9y 


The  curvature  /c  is  defined  as 
1 


(5) 


K  R ’ 

with  radius  E  shown  in  Fig.  3,  where  W  denotes  the  maximum 
displacement  of  the  multilayer  laminate  at  the  center. 

Assuming  each  layer  is  linear  elastic,  isotropic,  stresses  in  each 
layer  can  be  expressed  as  [24]: 

cr  =  Q(e-aATl ),  (6) 

where 

a  =  [crx,  cry,  rxy]T,  and  1  =  [1,1,0]T. 

and  a  is  the  coefficient  of  thermal  expansion  (CTE)  and  AT  is  the 
temperature  change.  The  above  equation  can  be  recast  into  the 
following  form: 


a  =  Q[£°  zk  a  ATI), 
where 


(7) 
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with  E  and  v  being  the  Young’s  modulus  and  the  Poisson’s  ratio. 
The  resultant  forces  and  moments  acting  on  all  the  layers  may 


be  calculated  as 

rh/2 

rh/2 

N  =  /  erdz, 

M  =  /  (rzdz 

(8) 

J -ft/2 

J-h/2 

with 

N  =  [NX,Ny,NXy]J, 

M  =  [MX,  My,  MXy]J. 

Substitution  of  (7)  into  (8)  yields: 
N  =  Ae°  -Bk-NqI, 

M  =  Be0  -  Die  -  M0 1, 
where 

"  rzi+ 1 

a=/2  Qdz’ 

k= i  ^ 


(9) 
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Qz2  dz. 


O';  ATE; 
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dz, 


r+1  aiATEi 

m°  =  E/  T^zdz- 

k=i  Jzi 

Note  that  for  functionally  graded  materials,  both  the  elastic  con¬ 
stant  matrix  Q  and  the  CTE  a  are  functions  of  the  coordinate  z. 
Furthermore,  if  the  temperature  distribution  is  not  uniform  across 
the  thickness,  AT  is  also  a  function  of  z. 

Consider  the  equilibrium  of  an  infinitesimal  element  in  the  lam¬ 
inate,  the  fact  that  the  total  forces  and  total  moment  acting  on  the 
element  must  add  up  to  zero  provides  the  following  equilibrium 
equations  [25]: 


dNx  3N; 


dNy  3N; 


^  +  ^  =  0,  + 


xy 


dx  dy 


dy 


dx 
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(11) 
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dx4 

=  -  (  q  +  Nx 


dx2 dy2  dy4 
92w0 


dx2 


dy2 


92w0 
'Xy3x3y 


(12) 


where  q  is  the  distributed  transverse  loading,  and  w0  represents  the 
z-displacement  of  the  midplane. 

With  various  boundary  conditions,  the  close  set  of  equations  will 
be  obtained  so  that  the  curvature,  displacement  and  stress  in  the 
various  layers  may  be  calculated.  The  several  boundary  conditions 
are  listed  as  follows: 

Free  edge:  For  a  free  edge,  the  bending  moment  and  the  shear  force 
must  vanish.  Thus, 


Mx(a,y)  =  0, 


dMx  t  n  dMxy  _  n 

~dx  ~W~ 


(13) 


Simply  supported  edge :  For  a  simply  supported  edge,  the  deflection 
and  the  moment  are  zero.  Thus, 


w0(a,y)  =  0,  Mx(a,y)  =  0 


(14) 


Rigidly  clamped  edge:  In  this  case,  the  deflection  and  the  rotation 
are  zero.  Thus, 


w0(a,y)  =  0, 


9w0(x,y) 


dx 


=  0,  e°(a,y)  =  0 


(15) 


2.2.  Global  fracture  criteria 


Assuming  a  crack/flaw  A  exists  in  a  thin  film  structure  as  shown 
in  Fig.  4,  thermal  loading  T  and/or  mechanical  loading  Twill  create 
a  corresponding  stress  field  o ,  where  U  represents  the  deformation 
field,  0  is  the  given  displacement  at  the  boundary.  G  refers  to  the 
energy  release  rate  at  the  crack  tips. 

Stress  increase  caused  by  sintering  and/or  the  stack  assembly 
processes  may  be  calculated  as 


cr=f(F,E,X), 


(16) 
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K  t  /* 


The  plate  curvature  due  to  temperature  differential  AT  can  be 
obtained  as  [24]: 


k  =  £2  AT, 


where 

£2  _  Nq(Bii  +  ^12)  -  M0(An  +^12) 
(Bn  +B12)2  -  (T>n  ±D12)(An  ±A12) 

with 


otjEj 
1  -  v 


dz, 


(25) 


where  E  and  X  represent  the  mechanical  properties  of  materials  and 
geometry,  respectively. 

The  energy  release  rate  at  the  crack  tip  is  denoted  as  G.  It  is  a 
function  of  both  the  mechanical  properties  E  and  the  geometry  X, 
such  that 

G  =  g{F,  E,  X).  (17) 

If  this  energy  release  rate  exceeds  its  critical  value,  Gc, 


G  =  g(F,  E,  X)  >  Gc, 


(18) 


M0 


1  -  v 


zdz. 


During  stack  assembly,  the  warped  cells  with  curvature  /c  are 
flattened  by  moment  M.  The  relationship  between  the  curvature  k 
and  the  flattening  moment  M  is  obtained  as 

k  =  T2M.  (26) 


The  corresponding  stress  increase  caused  by  cell  flattening  can 
now  be  calculated  as 


fracture  failure  occurs.  It  is  a  local  fracture  criterion  based  on  the 
stress  field  at  the  crack  tip  obtained  either  numerically,  analytically, 
or  experimentally. 

In  the  meantime,  the  thermal  and/or  mechanical  load  F  induces 
deformation  of  the  structure,  U,  as 

U  =  V{F,E,X).  (19) 

Combining  Eqs.  (18)  and  (19),  the  relationship  between  energy 
release  rate  and  deformation  may  be  obtained  as 

U  =  <I'(g~\C),E,X).  (20) 

If  the  relationship  given  by  Eq.  (20)  is  monotonous,  the  critical 
deformation  Uc  can  be  determined  given  the  critical  energy  release 
rate  Gc, 

Uc  =  ng-\Cc),E,X).  (21) 

Therefore,  fracture  failure  will  occur  if 

U  =  E,X)  >UC  =  g-’CGc ),  E,X).  (22) 

Hence,  the  fracture  criterion  is  based  on  the  global  variable,  i.e.,  dis¬ 
placement  field,  which  may  be  measured  easily  after  the  sintering 
processes  of  the  SOFC  PEN  structures. 


3.  Global  failure  criteria  for  SOFC  PEN 

3.1  Thermal  residual  stress  induced  deformation  of  SOFC  PEN 
structures 


During  the  sintering  process,  SOFC  PEN  structures  are  placed  in 
the  furnace  freely  without  any  constraints  on  the  boundary. 

Assuming  free  edge  condition  for  the  PEN  boundary  at  x  =  ±1/2 
and  y  =  ±1/2,  the  bending  moment  and  the  shear  force  must  vanish, 
which  leads  to  the  following  set  of  boundary  conditions  [25]  as 
boundary  condition  given  in  (13): 

Mx(±i,y)=°,  Mx(x,±P)=  0  (23) 


dMx  t  n  dMxy  —  n 
~dx~  +Z~dy~  ~  U 


(24) 


a  =  zQk 


(27) 


3.2.  Global  failure  criteria  for  SOFC  PEN 

Due  to  the  porous  nature  of  anode  and  cathode  in  the  PEN  struc¬ 
tures,  initial  flaws  and  crack  initiation  points  will  naturally  exist  on 
the  interfaces  of  anode/electrolyte/cathode  and  in  the  interior  of  the 
materials.  The  types  of  initial  flaws  considered  here  are  schemati¬ 
cally  illustrated  in  Fig.  5. 

Although  various  techniques  have  been  developed  for  sintering 
the  cell  materials,  certain  degrees  of  cell  warpage  will  remain  after 
sintering.  During  the  stack  assembly  process,  cells  are  flattened  in 
order  to  fit  into  the  stack.  Flattening  exerts  additional  stresses  to 
the  warped  cells.  Such  stress  may  fracture  the  cell  during  the  flat¬ 
tening  process  if  the  warpage  is  greater  than  some  tolerable  level. 
Therefore,  the  global  failure  criteria  developed  here  for  PEN  struc¬ 
ture  of  SOFCs  are  based  on  local  fracture  mechanics  formulations 
for  initial  crack  size  at  various  location  of  the  cell. 

For  most  anode-supported,  YSZ-based  SOFCs,  the  crack  types 
A,  B,  C  and  D  are  particularly  vulnerable  to  the  flattening  process. 
The  failure  criteria  below  provide  a  first-order  estimate  of  the  max¬ 
imum  allowable  warpage  and  maximum  allowable  curvature  to 
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Fig.  5.  Cracks  at  various  locations  in  a  cell,  where  A,  surface  crack  in  the  cathode;  B, 
surface  crack  in  the  anode;  C,  interfacial  crack  between  the  cathode  and  electrolyte; 
D,  interfacial  crack  between  the  anode  and  the  electrolyte;  E,  blister  crack  on  the 
anode/electrolyte  interface;  F,  tunneling  crack  in  the  electrolyte 
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avoid  fracture  during  the  stack  assembly  process.  In  the  follow¬ 
ing  context,  a  is  the  crack  length  of  each  crack  type,  Ea(e,c>  fa(e,c> 
aa(e,c)  and  fra(e,c)  represent  the  Young’s  modulus,  Poisson’s  ratio,  CTE, 
and  thickness  of  the  anode,  electrolyte,  and  cathode,  respectively. 
In  addition,  subscript  ce  and  ae  stand  for  the  combined  properties 
of  cathode/electrolyte  composite  and  anode/electrolyte  composite, 
respectively. 


curvature  /c  is  flattened  by  moment  M,  the  relationship  between  /c 
and  M  can  be  obtained  in  Eq.  (23)  as 


M(An  +AU) 


(Bn  +B12)  -  (£>n  +Du)(An  +  A12) 


(34) 


Assuming  free  boundary  conditions  at  the  edge,  the  average  stress 
in  the  cathode  may  be  obtained  as 


3.3.  Crack  A 


a  = 


h  - 

~^2 


a 


Eck, 


(35) 


First,  we  consider  the  critical  curvature  the  tri-layer  PEN  struc¬ 
ture  can  withstand  during  cooling  from  the  sintering  temperature, 
i.e.,  stress-free  temperature  to  room  temperature  given  the  initial 
flaw  A. 

Based  on  elasticity  theories,  the  average  biaxial  stress  in  the 
cathode  and  electrolyte  for  temperature  differential  AT  can  be 
approximated  as  [22]: 


where  h  represents  the  thickness  of  the  PEN  structure.  Substitution 
of  Eq.  (35)  into  (30)  leads  to 

G=(^)2£cHV.  (36) 

Superimpose  the  energy  release  rates  from  sintering  and  flat¬ 
tening,  i.e.,  combining  Eqs.  (33)  and  (36),  one  obtains 


a  =  — —AT Ace, 
1  -  vc 

where 


(28) 


A  a  =  ac  -  Oiae, 


Oiae  — 


OtaEaha  +  OieE^he 
Oiah-a  T-  OCehe 


Noted  that  AT  is  the  temperature  difference  from  the  stress-free 
temperature  to  room  temperature.  Combining  Eqs.  (25)  and  (28) 
yields 


k  =  Q 


1  -  Vc 
EcAa 


cr, 


(29) 


For  stress  a  in  this  crack  configuration,  the  energy  release  rate  may 
be  approximated  by  the  following  relations  [21  ]: 


G  = 


K2 
£c  ’ 


(30) 


where  I<,  the  stress  intensity  factor,  is  approximated  for  this  type  of 
crack  by 


K  =  crH 

H.1.,2,5V5i(l-^)',“(l+AiL). 


Here  a  is  the  initial  crack  depth,  A  is  a  fitting  parameter  which  is 
independent  of  a/h3  [21  ],  and  s  is  the  root  of  the  following  equation 
[26]: 

cos(s7r)-2^^(l-s)2  +  ^^=0.  (31) 

1  -  P 2  1  - 

where  j3\  and  are  the  two  Dundurs  parameters  following  the 
work  of  Dundurs  [27].  The  relationship  of  the  critical  energy  release 
rate  Gc  and  critical  curvature  kc  is  obtained  as 


Cc  = 


1  A  a 

£2  (1  -  vc) 


EcH2k2, 


(32) 


Cc  = 


'  1  Aa 
.  £2  1  -  vc  + 


(t2) 


EcHzk ; 


2,^2 


(37) 


For  given  critical  energy  release  rate  Gc,  the  critical  curvature  /cc, 
is  then 


'  1  Aa  fh-a\~  1  1  /g7 

.32  1  -  vc  +  V  2  )\  H\Tc' 


(38) 


3.4.  Crack  D 


Considering  the  edge  crack,  Crack  D,  located  in  the  interface 
of  the  anode  and  electrolyte  in  the  PEN  structures.  Again,  crack 
tip  field  generated  by  cooling  is  considered  first.  The  cathode  and 
electrolyte  in  the  PEN  structure  is  considered  as  bi-layer  compos¬ 
ite  plate.  For  bi-materials  laminate,  the  equivalent  elastic  modulus, 
Poisson’s  ratio,  and  coefficient  of  thermal  expansion,  can  be  given 
as  [28]: 


heEe  +  hcEc 
Ece 


fee 


EqVq  H-  YiqVq 
Ece 


a^heE^  -\-  achcEc 

heEe  +  hcEc 


(39) 

(40) 

(41) 


The  energy  release  rate  may  be  obtained  under  the  temperature 
differential  as  [16]: 


G=  — F1(Aa)2(AT)2, 


(42) 


with 


C 2d2  d2 

E\  =  1  1  +  -jr  +  ^2C3d\d2d3. 


The  parameters  in  the  above  equation  are  defined  in  Appendix 
A  as  Ref.  [26]. 

From  Eqs.  (25)  and  (42),  one  may  obtain  the  relationship 
between  the  energy  release  rate  and  the  curvature  as 


and 


£2  1  -  vc  l~Gt c 
Kc=  H  Aa  Y  Tc 


(33) 


For  given  critical  energy  release  rate  Gc,  the  critical  curvature 
can  now  be  determined. 

Next,  the  energy  release  rate  at  the  crack  tip  generated  by  the 
flattening  process  during  stack  assembly  is  superimposed  on  that 
resulted  from  the  cooling  process.  Assuming  a  warped  cell  with 


^  1  4hCeF\(Aa)2  2 

322  Cce 
therefore, 


k  =  Q 


Cce 

4hceF(Aa)2 


1/2 

Vc. 


(43) 


(44) 


Next,  crack  tip  field  resulted  by  cell  flattening  during  stack 
assembly  is  derived.  For  a  cell  with  curvature  /c,  the  relationship 


W.N.  Liu  et  al.  /  Journal  of  Power  Sources  192  (2009)  486-493 


491 


between  the  crack  tip  energy  release  rate  and  the  curvature  and  is 
expressed  as 


Ccef2  K2 

16  h3ce  &  ’ 


(45) 


which  leads  to, 


k  =  £2 


16  h3eG 
CceF2 


1/2 


(46) 


where 


_ +^12 _ 

(Bn  +B12)2  -  (£>n  +^12X^11  +  A\2) 


F2 


cl  +  cl 

r  1 


+  12C2C3d3. 


Again,  combining  the  effects  of  cooling  and  cell  flattening,  the 
relationship  of  the  energy  release  rate  and  the  curvature  can  be 
obtained  as 


1  ( EsLl]  1/2 , 1  / jhceMAaf 

£2  \  16hce  /  Cce 


(47) 


Using  similar  derivation  process,  for  each  crack  type  illustrated 
in  Fig.  5,  one  may  derive  the  relationship  between  the  critical  energy 
release  rate  and  critical  curvature  of  the  warped  PEN  structure  as 


kc  =  Y 


(48) 


where  Y  is  dimensionless  constant  depending  on  the  cell  geometry, 
material  constants  as  well  as  the  size  and  location  of  each  crack  type. 
For  Crack  A, 


Y  = 


'  1  Aa 
_  £2  1  -  vc  + 


heEe 


and  for  Crack  D, 


1  f  cCeE2  \  1  /  4hceF\{Aa) 

J2\16h3ce)  c^e 


1/21 


-1 


\J  heEe- 


(49) 


(50) 


Similarly,  the  relationship  between  the  critical  energy  release 
rate  and  maximum  deformation  of  the  warped  PEN  structure  may 
be  obtained  as 


where  W  is  the  maximum  allowable  warpage,  L  is  the  nominal  cell 
size,  Y,  like  Y,  is  dimensionless  constant,  too,  depending  on  the  cell 
geometry,  material  constants  as  well  as  the  size  and  location  of  each 
crack  type. 


4.  Illustrative  examples 

As  an  example,  consider  a  cell  with  geometry  and  material 
parameters  listed  in  Table  1.  For  this  example,  the  normal  stresses 
are  tensile  in  the  cathode,  and  compressive  in  the  electrolyte  and 


Table  1 

Material  properties  and  geometric  parameters  of  a  cell. 


Table  2 

Values  of  the  factor  Y  ( x  10-3 ). 

Crack  A  Crack  C  Crack  D 


a  =  0.01/i3  a  =  0.05/i3  a  =  0.1h3  a  =  0.2/i3 


0  5  10  15  20 

Gc  (J/m2) 

Fig.  6.  kc  vs.  Gc  for  the  cell  parameters  given  in  Table  1. 


anode.  Therefore,  crack  types  A,  C  and  D  may  propagate  during  cell 
flattening.  The  values  of  Y  for  these  crack  types  are  listed  in  Table  2. 
For  L  =  10  cm,  the  relationship  between  the  curvature  and  the  energy 
release  rate  is  depicted  in  Fig.  6  for  various  crack  types;  similarly, 
Fig.  7  illustrates  the  relationship  between  the  maximum  allowable 
warpage  W  and  the  fracture  toughness  Gc  for  these  cracks. 

It  is  clear  that  the  critical  curvature  and  maximum  allowable 
deformation  for  the  PEN  structure  with  the  surface  crack  on  the 
cathode  will  be  lowered  with  the  increment  of  the  crack  length 
with  given  critical  energy  release  rate.  Compared  with  the  interfa¬ 
cial  edge  crack  between  the  anode  and  the  electrolyte,  the  critical 
curvature  and  maximum  allowable  deformation  for  the  interfacial 
edge  crack  between  the  cathode  and  electrolyte  is  much  less.  The 
interfacial  edge  crack  between  the  cathode  and  electrolyte  and  the 
channeling  crack  in  the  cathode  will  be  easier  to  result  in  fracture 
failure  of  the  PEN  structure  after  the  sintering  process  and  during 
the  cell-stack  assembly  process.  The  crack  and  flaw  size  in  the  PEN 
structure  may  be  controlled  by  the  sintering  processes  parameter 


Fig.  7.  W  vs.  Gc  for  the  cell  parameters  given  in  Table  1. 
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Fig.  8.  Relationships  between  critical  curvature  and  displacement  and  critical 
energy  release  rate. 
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Appendix  A 


Cce  — 


Yce  +  1 
fLce 


Yce  —  3  —  4vce, 


Y\  =  3  -  4vi , 


hce  —  fl2  +  fl3 , 


_  l^ceiYl  +  2 )  ~  l^S^Yce  +  1 ) 

A^ceC/l  +  1)  +  1^3  (Yce  +  1) 


and  the  component  materials.  The  process  parameter  and  compo¬ 
nent  materials  can  be  optimized  with  combination  of  the  fracture 
failure  analysis  developed  here  in  order  to  gain  the  reliable  PEN 
structure. 

For  each  crack  existing  in  the  PEN  structure  of  SOFCs,  for  a  given 
interfacial  strength  (i.e.,  interfacial  fracture  toughness  Gc),  the  cor¬ 
responding  critical  curvature  kc  can  be  determined.  This  critical 
curvature  and  maximum  allowable  deformation  of  the  PEN  struc¬ 
ture  will  help  to  judge  if  PEN  structure  survives  during  the  cell-stack 
assembly  process  pursuing  after  the  sintering  process.  For  example, 
as  shown  in  Fig.  8  for  Crack  C,  if  the  curvature  and/or  the  maximum 
deformation  of  the  PEN  structure  induced  by  the  sintering  process 
are  larger  than  the  critical  value  obtained  here,  fracture  failure  will 
occur  during  cell-stack  assembly  (flattening).  Therefore,  curvature 
and/or  warpage  of  the  PEN  structure  may  be  measured  after  the 
sintering  process,  the  PEN  structure  has  to  be  abandoned  if  their 
values  of  the  curvature  and/or  warpage  are  larger  than  the  critical 
value  obtained  here. 


$  = 


\  +  A 
1  -  A  ’ 


5  = 


l+2gyy  +  gq* 2 
2*7(1  +  %rj) 


Ao  =  — E 
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5.  Conclusions 

In  this  paper,  a  suite  of  global  failure  criteria  for  the  SOFC  PEN 
structure  is  derived  based  on  fracture  mechanics  formulations  of 
local  crack  tip  field  for  cracks  at  different  locations  with  different 
initial  sizes.  These  failure  criteria  consider  the  stresses  generated 
from  the  initial  PEN  cooling  from  sintering  as  well  as  the  mechanical 
flattening  process  during  stack  assembly.  The  purposes  of  develop¬ 
ing  such  global  failure  criteria  are  to  aid  the  initial  design,  material 
selection  and  optimization  of  SOFC  fabrication  process.  The  critical 
energy  release  rate  is  related  to  critical  curvature  and  maximum 
deformation  of  the  warped  cells.  Since  the  curvature  and  displace¬ 
ment  of  the  warped  cells  can  be  measured  prior  to  stack  assembly, 
engineers  can  use  the  global  failure  criteria  developed  here  to  pre¬ 
dict  possible  failure  of  a  cell  structure  for  various  initial  crack  types. 
Such  global  failure  criteria  can  also  be  used  for  designing  the  sinter¬ 
ing  processes  such  that  allowable  warpage  and  curvature  level  of 
the  PEN  can  be  achieved  for  its  survival  in  the  subsequent  assembly 
process. 
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